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Structure of cold dense stellar matter at subnuclear densities is investigated with quan¬ 
tum molecular dynamics (QMD). We succeeded in showing that the phases with slab-like 
and rod-like nuclei etc. can be formed dynamically without any assumptions on the nuclear 
shape. Our result suggests the existence of these kinds of phases in neutron star crusts. 

§1. Introduction 

Collapse driven supernovae (SNe) and the following formation of neutron stars 
(NSs) are the most dramatic processes during the stellar evolution. These objects 
provide not only astrophysically significant phenomena but also interesting material 
phases inside them; both are strongly connected with each other. At subnuclear 
densities in NS crusts and SN cores, nuclei are expected to have exotic structures 
such as rod-like and slab-like shape etc., which me called nuclear “pasta”. HP This 
prediction has been obtained by previous studiesIHP’Sl based on free energy calcula¬ 
tions assuming specific nuclear shapes. These works clarify that the nuclear shape 
is determined by a subtle balance between nuclear surface and Coulomb energies. 

However, the formation process of the “pasta” phases has not been di^iissed 
except for some limited cases which are based on perturbative approaches. HP It 
is important to adopt microscopic and dynamical approach which allows arbitrary 
nuclear structures in order to understand the formation of the non-spherical nuclei. 
We have started studying dense matter by the QMD HP simulation whose final goal 
is to obtain a realistic picture about NS crusts and SN cores. As the first step of our 
quest, we will answer the following question in this study: “Can the ’pasta’ phases 
be formed dynamically?” 


§2. Simulations and results 

In the formation process of the “pasta” phases, density fluctuations would be 
crucial. Methods of treating fermion many body systems based on molecular dynam¬ 
ics are suitable for incorporating fluctuations of particle distribution. We choose the 
QMD among a lot of versions of MD for fermions based on trade-off between the 
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Fig. 1. Nucleon distibutions of cold symmetric nuclear matter at various densities; (a) sphere phase, 
O.lpo; (b) cylinder phase, 0.25po; (c) slab phase, 0.4po; (d) cylindrical hole phase, 0.5po and 
(e) spherical hole phase, 0.6po. The red particles represent protons and green ones represent 
neutrons. Temperatures are 0(0.1) MeV. 


calculational amounts and accuracies. It is noted that we watch macroscopic struc¬ 
ture, for which the exchange effect is not so important. Therefore, it is justified to 
use the QMD which is less elaborate in treating the exchange effect. 

We have performed QMD simulations of (n,p,e) system with proton fraction 
X = 0.3 and 0.5 for various nucleon densities p. We set a cubic box which is imposed 
periodic boundary conditions. Total nucleon numbers are 1372 and 2048. The 
Coulomb interaction is calculated by the Ewald method and relativistic electrons 
are treated as a uniform background. For lie nucleon interaction, we use a model 
Hamiltonian developed by Maruyama et al.0) 

We first prepare uniform hot nucleon gas at ksT ~ 20 MeV as an initial condi¬ 
tion, and then cool it down slowly for 0(10^ — 10^) fm/c until the temperature gets 
0(0.1) MeV. Note that any artificial fluctuations are not given during the simula¬ 
tion. Shown in Fig. 1 is the resultant nucleon distributions of cold symmetric matter 
for p = 0.25 — O.Gpo ( Po • the normal nuclear density ). We can see from Fig. 1 
that the phases with rod-like and slab-like nuclei, cylindrical and spherical bubbles, 
in addition to the one with spherical nuclei, are reproduced. Also for x = 0.3, these 
phases are reproduced at about 0.05/9o lower than the symmetric matter case. 

The density dependence of the nuclear shape is quite sensitive and various inter¬ 
mediate phases are observed in density regions between the densities corresponding 
to the phases with the simple shapes as shown in Fig. 1. Whole phase diagrams of 
cold matter and properties of the nuclear shape changes will be discussed elsewhere.™ 
It is concluded that our results suggest the existence of the “pasta” phases in NS 
crusts which cool down keeping the local thermal equilibrium after the proto-NS is 
formed, and their cooling time scale is much longer than that of our simulations. 
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